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I. INTRODUCTION 



The dramatic discovery of Charmonia, the J /if) and its excited states, marked the 
beginning of a new era of particle physics. Till now the Charmonium physics remains to 
be one of the most exciting areas of high energy physics. As the " hydrogen-like atoms" of 
strong interaction, the Charmonia could be investigated partly by virtue of the perturba- 
tive QCD (pQCD) in account of the large charm quark mass, which makes the study on a 
relatively solid ground, as well the study may give clues of the nature of non-perturbative 
QCD. 

Although the first Charmonium state, J/if), was observed more than twenty years 
ago, the study of the Charmonium states is still far from satisfactory. Except for the 
J /if) itself, the knowledge of the other Charmonia is very limited. We do not even have 
a complete cc mass spectrum below the DD threshold JTJj, that is, the existence of the 
S-wave spin singlet i]' c and the P-wave spin singlet h c {}Pi) is still based on very weak 
experimental signals. To confirm the existing findings and give out more precise values of 
the mass, width, and other parameters of these two resonances are now a pressing task 
in experiment. 

The r]' c was first observed in the Crystal Ball experiment in the inclusive photon energy 
spectrum from ip' decays at 3594 MeV || , until now the signal was not observed by other 
experiments due to the low energy of the radiative photon and the relatively poor photon 
detection ability of other detectors comparing with that of the Crystal Ball's The 
h c i}Pi) state was first observed at 3526.14 MeV in the proton antiproton annihilation 
experiment by E760 group performed at Fermilab 0, but with the statistical significance 
of the signal slightly more than three standard deviations, and also no other experiments 
definitely confirm the existence by now (the E705's report [Q was doubted by Barnes, 
Browder, and Tuan 

Currently, the experiments suitable for Charmonium studies are: the BES detector 
running at the BEPC e + e~ collider, the pp annihilation experiment represented by the 
E835 experiment at Fermilab, and the scarce studies of the two-photon process in high 
energy e + e~ colliders like at LEP and CESR. 

Due to the restriction of the quantum number at e + e~ colliders, only the vector-like 
Charmonium states like J /if) and if>' can be produced directly at lowest order, whereas 
the other Charmonium states, like Xcj, T] c , and /i c ( 1 Pi), can only be produced via either 
higher order processes or through the J/if)(if)') electromagnetic and/or hadronic decays. 
For instance, the v/ c may be produced via if)' — > 7 + r)' c and the h c (}Pi) state via if)' — > 
7T° + h c (}Pi). Although BES detector |T0| has collected the largest if)' data sample in the 
world, due to the limited energy resolution of the Electromagnetic Calorimeter and the 
rather small production rates of r/ c and h c {^P\) in if)' decays, the search of either v/ c or 
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/i c ( 1 Pi) did not give significant results. 

As for proton-antiproton annihilation experiments, although they can produce Char- 
monium states of various quantum numbers and can be used to determine the resonance 
parameters of the Charmonium produced, the study of Charmonium is limited by the 
detection of the electromagnetic final states and the low production rate. The E760 and 
its succeeding version E835 did a very good job in measuring the resonance parameters 
of the Xci, Xc2 and some other Charmonium states, but the study on the h c (}P\) state 
is still insufficient and the existence of the rj' c is not confirmed. Of course, the E835 will 
continue this work and look further with more data in the near future. 



The HERA-P an experiment presently set up at DESY, which uses the HERA 
920 GeV proton beam incident on various nuclear targets, is focused on the measurement 
of CP- violation in the BB system via mainly the final states containing J/ip. The trigger 
system is designed to recognize events with J/ip — > £ + £~ (£ = eor/i). Furthermore, 
the detector also is designed for precise measurement of photons with its Electromagnetic 
Calorimeter (ECAL), which makes the study of Charmonia very possible through detecting 
the final states of the Charmonium decays containing J/ip and neutral particles like 7 or 



7T°. 



The paper is organized as follows. In following section we present the formalism for 
r)' c and /i c ( 1 Pi) production in a general framework in fixed-target experiment. In section 
III the obtained formalism is applied to HERA-P situation numerically; the direct and 
indirect production rates of these two states are evaluated. In section IV, we give a 
rough estimation of the signals and backgrounds in searching for these two states for 
experimentalists' reference. In the last section some discussions and conclusions are made. 



II. rf e AND h^Px) PRODUCTION 



For rf c production, to leading order in a s and v 2 , the relative velocity of heavy quarks 
inside the bound state, it is a two to one process as shown in Figure 1. The parton level 
cross section can be easily calculated or just obtained from the corresponding r\ c producing 
process with the non-perturbative sector replaced. It is 

" 1 = 972^ <0| ^ (l5o)|0> ^ (1 -" ) - (1) 

Here, a s is the strong coupling constant; < 0\O^ c ( l S )\0 > is the NRQCD Color-Singlet 
non-perturbative matrix element, which can be related to |_R^(0)|, the radial wave func- 
tion at the origin of the bound state, by < OlO^ (} S )\0 >= ^|-R(0)| 2 ; and z = M^Js, 
where s denotes the c.m.s. energy in partonic system. 
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As for the h c (}Pi) production, the situation is somewhat different from that of r)' c . 
Of the latter, at leading order in a s and v 2 there is only one possible channel giving the 
contribution, but of the former, there are several to the same order of accuracy. To be 
more clearly, according to the BBL theory for Quarkonium production and decays ||12|| , 
the Fock states of Quarkonium are ordered in v, i.e., 

I^CPO) =0(1) caM iy ]) + G(v) cc[ 1 S^)g) + 0(v) cc^f ] g) + ■ ■ ■ . (2) 



Because for P-wave states the leading non-vanishing wave functions are the derivative of 
the wave functions at the origin, or in other words that the P-wave states are produced 
via the NRQCD dimension 8 operators or higher, the NRQCD scaling rules fl3| tell us 
that for /i c ( 1 Pi) production the non-perturbative matrix elements stemming from the first 
two terms in Eq.(fJ) are of the same order in v 2 . Based on this argument the leading order 
Color-Singlet and -Octet processes of the h c i}P\) production are shown in Figure 2. 

As depicted in Figure 2(a), the Color-Octet process is also a two to one process. The 
cross section of the partonic scattering process can be straightforwardly obtained, 

°2 = yI^^ < O|0fcf£ o )|O > - z) , (3) 

where < 0|0 8 c ( 1 5'o)|0 > is the Color-Octet nonperturbative matrix element. 

Of the Color-Singlet processes, Figure 2 (b)-(d), the two gluon fusion channel of (b) 
may survive only with at least an additional gluon in the final states from the Landau- 
Yang theorem, as shown in the Figure; the others are not restricted by this law, however 
ruled out by the properties of charge-conjugation of the processes. The reason for this 
is that heavy-quark-loop factor (including the projector for the quarkonium state) is odd 
under charge conjugation. That is, the C-odd h c state can not decay through two vector 
currents (C-even), and the direct calculation really shows they give no contributions. The 
cross section of Figure 2(b) reads as 

H9 + 9^ h^P^]) _ rc 2 a 3 < Offi ? P^O > / 2 ^t 2 u 2 + stu{t + u) + 2s 2 (t 2 + tu + u 2 



dt 108(2m c )s 2 [ (s + i) 2 (s + u) 2 (t + u) 2 

40 

+ ; s (12s 6 t + UsH 2 + 72s 4 t 3 + 72s V 

+ Us 2 i 5 + 12s£ 6 + 12s 6 m + 58sHu + U9s 4 i 2 u + 179s 3 t 3 u 

+ U0s 2 i 4 u + 56si 5 u + 12t e u + A6s 5 u 2 + I57s 4 iu 2 + 2A6sH 2 u 2 

+ 231s 2 t 3 u 2 + U2si 4 u 2 + Ut 5 u 2 + 78s 4 m 3 + 198s 3 iu 3 + 240s 2 Pn 3 

+ 178si 3 u 3 + 72t 4 u 3 + 79s 3 u 4 + 158s 2 M 4 + U9si 2 u 4 + 72t 3 u 4 

+ 47s 2 m 5 + 61siu 5 + Ut 2 u 5 + 12su 6 + 12tu 6 ) } . (4) 

Here in the above, the s = (P1+P2) 2 , t = (pi— P3) 2 , and u = (P2—P3) 2 are ordinary Mandel- 
stam variables; the universal non-perturbative matrix element < 0|(9j' c ( 1 Pi)|0 > related 
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to the derivative of the radial wave function at original of h c ( lp i) h Y < 1 (9^ c (^i) |0 >= 

£K c (o)l 2 - 

Except for the direct production of these two states given in above, another main source 
of their production is of the electromagnetic or hadronic decays of the ip' in accompanying 
with one 7 or 7T°. The dominant partonic interaction processes of the ip' production in 
pN collision at HERA- 5 energy are drawn as Figure 3. 

The expression for gluon-gluon fusion processes, the Figure 3(a) and (b), can be written 

as 

^4(9 + 9 -> = 

{< 0\OfCS )\0 > + ± < 0\Ot'eP )\0 > < 0|^'( 3 P 2 )|0 >} z6(l - z) 

+ 20n2a i « oi<(^)io > , 2 r~f +2z i ogz + 1- f +2 :! ogz ]e(i - Z ) . ( 5) 



81(2m c ) 5V 1 1 v 1/1 [ {z-lf (z + lf 



The expression for process of Figure 3(c) is quite simple, it is 

lQ7r 3 a 2 s 
27(2m c ) 5 



^{q + q - n 3 S[ 8) ]) = n ^f a t < 0\Oi'CS )\0 > z5(l - z) . (6) 



Here, although the Octet processes are suppressed in v 2 , they get compensation from 
the enhancement of l/a s relative to the Color-Singlet process. So, it is proper to include 
them in the ip' production rate estimation. 



III. NUMERICAL ESTIMATION FOR rf c AND h^Px) PRODUCTION AT 

HERA-B 

In the above section we have calculated the necessary partonic cross sections at lead- 
ing order in v 2 or/and a s for r)' c and /i c ( 1 P 1 ) production in the proton-nucleon collision. 
According to the general factorization theorem the experimental cross sections can be 
obtained by convoluting the subprocess with the parton distribution functions in the 
nucleons. i.e., 

a(A + B^C + X) = J2j G a (x a )G b (x b )a(a + b^C + Y)dx a dx b , (7) 

where the sum runs over all the possible initial interacting partons which involve in the 
interaction; the A and B represent nucleons; C represents the Charmonium; X and Y 
are the remnants of the inclusive processes; G a (x a ) and G b (x b ) are the parton distribu- 
tion functions of the colliding nucleons A and B with momentum fractions x a and x b , 
respectively. 
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In doing the numerical estimation the following inputs are taken 



a s {2m c ) = 0.253, M v , c = 3.6 GeV, M hc{ i Pl) = 3.5 GeV, m c = 1.5 GeV, 

< 0\O^CS )\0 >= 0.98 x 10~ 2 GeV 5 [14], < Offi C Pi)\0 >= 0.32 GeV 5 |5| 



< 0\ofCS )\0 > +^ < 0|Of'( 3 P )|0 >= 0.56 x 10~ 2 GeV 3 p 
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m 



< 0|Cf ( 3 5i)|0 >= 0.44 GeV 3 0, < 0|Of ( 3 Si)|0 >= 6.2 x 10~ 3 GeV 3 0, 

< 0|Of (^q)^ >= 0.20 GeV 3 p 



and the CTEQ 3M package for parton distributions is employed with the factorization 
scale chosen to be equal to the NRQCD scale /i = 2m c . In making use of the present 
fitted matrix elements given in above, the spin symmetry relation < 0\O% ( 3 Pj) |0 >= 
(2J+ 1) < 0|Of ( 3 P )|0 > has been applied. 

With 920 GeV incident proton we find the magnitude of the cross sections given in 
the preceding section are 



d! = 1076.1 nb/n, a 2 = 98.9 nb/n, a 3 = 54.8 nb/n, a 4 = 79.0 nb/n, a 5 = 5.2 nb/n . (9) 



Here, the nb/n means nb/nucleon for shorthand. The ip' production cross section 
(84.2 nb/n) agrees well with the experimental measurement of (75 ± 5 ± 22) nb/n by 
E789 ||19|| , indicating the reliability of the other calculations in this paper. However, 
quarkonium production rates are often sensitive to the choice of m c and the parton dis- 
tributions. To see the effect of the former, we assume the difference between calculated 
and measured i/j' production cross sections is a pure effect of m c , to cover the error of the 
measured value, m c should vary from 1.45 to 1.65 GeV. By changing m c from 1.5 to 1.45 
and 1.65 GeV in all other cross section calculations, the relative uncertainties of the as 
are shown below. As for the latter, we simply take another parton distribution functions, 
the GRV [pOfl , the deviations of the as are also listed below. 

Aax =l*°i + 28.5 % , Aa 2 =1^.3 + 2g 5 % ^ = +32.6 + 28>g % ^ 

Aa 4 =t 2 i 3 7 + 29.2 % , Aa 5 =1^4 - 5.8 % . (10) 



Here, the first deviations come from the the change of m c ("+" for m c = 1.45 GeV and "-" 
for m c = 1.65 GeV); the second corresponds to the choice of a different parton distribution 
code (GRV results relative to the CTEQ ones); We can see from the above results that 
the deviations of the cross sections relative to different parton distributions agree within 
30%, and the charm quark mass uncertainty changes cross sections around 50%. 

Due to the projected high interaction rate, 40 MHz, the results in Eq. (|9[) means 
that in a running time of 10 7 s at HERA-S using the Cu target, for example, the directly 
produced r\' c and /i c ( 1 Pi) events number would be about 3.3 x 10 10 and 4.7 x 10 9 . The ip' 
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events number would be about 2.6 x 10 9 , which is three orders higher than the present ij)' 
date sample collected at e + e~ colliders. 

Theoretical estimation of the branching fractions of the /i c ( 1 Pi) production in ip' decays 
are about 3 from Refs. [f2~l"| -|2"3|l , and the i]' c rate are about 10~ 4 3 from the naive 

estimation of the Ml transition in non-relativistic limit [23,25] . Therefore, the indirectly 
produced /t c ( 1 Pi) and r( c would be of the order 10 4 ~ 6 and 10 5 ~ 6 correspondingly. 

The indirect production of Charmonium in B decays has been estimated in Ref. ||26|| , 
the production rates of h c i}P\) and rf c (assuming the same as that of r/ c ) are of the order 
of 10~ 3 . Using the bb production cross section of 12 nb/n, the produced hJ^Px) and r]' c 
events are of the order of 10 5 ~ 6 in 10 7 s of the HERA-P running time, which is the same 
order as via ip' decays. 



IV. SEARCHING STRATEGY 

As mentioned in the introduction part of this paper, the interested topologies of de- 
tecting these two states at HERA-B are ■yJ/ip and ir°J/i/; for T]' c and h^Pi) respectively, 
where J/ %fj decays into lepton pairs and 7r° decays to two photons. Because of the charge- 
conjugation invariance, the decay modes rj' c — > Tr°J/tp and h c — > "fJ/ip are ruled out. 

The /i c ( 1 Pi) state was observed decaying to 7i J/ip with branching ratio ~ 10~ 3 0, 



which is of the same order of magnitude as the theoretical expectation |2I|, and the rj' c 



decaying to jJ/ift is expected with a width of the order ~ Oil keV) p7| . Considering that 



the theoretical estimation of the decay width of the rj c is about 5 MeV, it has a branching 
ratio of ~ O(10~ A ) in jJ/ip decay mode. Using the numbers listed above, TABLE | lists 
the estimation of produced events for hd^Pi) and r)' c in all the production mechanisms, 
taken into account the branching ratios of J/ip leptonic decays and tt° — > 77. 

From the table, we can see the produced events of interested topologies from indirect 
productions are too low (of the order of 10 to 100) to produce meaningful signals for 
observing the two states. But instead, the direct productions of these two states are 
rather large, of the order of 4 ~ 6 x 10 5 . As we know the geometric acceptance of HERA- 
B detector is large and its trigger is optimized for J/ ip events, we do expect high efficiency 
of detecting these two final states. Suppose the overall efficiency of detecting these two 
final states is around 10%, one expects 4 ~ 6 x 10 4 reconstructed events each channel, 
which are large numbers compared to those channels for observing CP violation (in the 
same running time, the reconstructed events of J/ipK s is estimated to be around 1400!). 

The main background channel for 7]' c observation is Xc2 lJ /ip, which has the same 
final states but much larger cross section and very near the expected i]' c mass. Using the 
measured cross section of \c2 by E771 |28[[ , the number of reconstructed Xc2 events is 
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estimated to be around 10 8 (the combinational background at Xc2 mass region is about 



the same size as Xc2 events as shown in Ref. PBfl). The significance of the observed r)' c 
depends strongly on the mass resolution of jJ/ip system and the mass difference between 
Xc2 and r)' c . Theoratical estimations of the rf c mass ranges from 3589 to 3631 MeV f29fl , 
and only experimental hint is at mass of (3594 ±5) MeV. For a 3.6 GeV mass rj' c , if the 
mass resolution is around 10 MeV or less, rj c will produce a long tail at high mass side of 
Xc2, and at mass higher than 3.6 GeV, the events is almost free from Xc2 background. If 
the mass resolution reaches 15 MeV or even larger, it will be hard to distinguish r]' c from 
Xc2- A larger m n i obviously will increase the possibility of resolving T]' c from the Xc2 tail, 
while a low mass rf c will more depend on the mass resolution. 

For /z c ( 1 P 1 ), the main background is from the Tr 7r J/ip produced by ip' decays. Com- 
pared with that in r)' c case, here the /i c ( 1 Pi) is at the phase space limit of TT J/ip system 
produced from if)' decays and the cross section of the latter is smaller than Xc2 by at least 
a factor of 3%. Furthermore, there is no other nearby resonance decays to the same final 
states. All these make the observation of h c (^P\) easier than rj' c . 

At the point of data analysis, for r]' c , instead of using the invariant mass of J/ip and 
the detected 7, using the mass difference between the £ + £~'y system and the £ + £~ system 
would be better in finding the signal, since the latter can compensate some of the effects 
due to energy losses of radiation and bremsstrahlung of the lepton tracks. 

In searching for the /i c ( 1 Pi) state, the reconstruction of the 7r° is also important for 
the event selection, and it is also a very good constraint to lower the background level 
greatly. As in the ^J/ip case, the mass difference method will be helpful to this channel 
as well. 

Finally, to check the results, the sideband method maybe useful. In both cases the 
J/ip mass sidebands, and in /i c ( 1 P 1 ) searching the tt° mass sidebands will tell us the shape 
of the background. The absence of the same peak in the mass spectrum of sidebands 
events will be a demonstration that the selection is reasonable. 

It is important to note that all above discussions are based on a sample of 10 7 s running 
time. With more statistics, instead of reconstructing photon from ECAL, one can detect 
converted photon to reconstruct rj' c and ^i c ( 1 Pi), as has been indicated by E771 [28|. In 



this case, the mass resolution will be significantly improved (5.2 ± 2.0 MeV for 7J/V 
system in E771 experiment), rj c will be resolved from Xc2 even it has a small mass. 

V. DISCUSSIONS AND CONCLUSIONS 

In this paper we have discussed the physics potential of HERA-P in detecting the 
r)' c and /i c ( 1 Pi). Our numerical results reveal that there are about 10 10 and 10 9 of r)' c 
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and /i c ( 1 Pi) events would be produced at HERA-P in 10 7 s of running time. A rough 
estimation shows that /i c ( 1 Pi) will be observed clearly in its it°J/iIj decay mode, and 
i]' c will be observed as a shoulder at high mass side of Xc2 m iJ/ip channel if the mass 
resolution is not too large. 

The searching strategies of these two states at HERA-P are given. The major back- 
grounds in the detection and the possible detecting measures are also discussed. 

It should be mentioned that the theoretical basement of our calculation in this paper, 
the NRQCD factorization, may not work well in the inclusive quarkonium production at 
full phase space, that is at small pr (pt not much greater than Aq CD ) region, which would 
cast some shadow on the validity of the results of the inclusive fix-target calculations. 
However, at least from our calculation on ijj' production, which agrees with the experiment 
value quite well, we are convinced to a certain degree of our other calculations in this 
paper. 

Last, it should be noticed that although the study proceeded in the paper is just an 
order estimation because either input parameters, like the color-octet matrix elements, 
are more or less accurate just to an order, or the evaluation is based only on the first 
order calculation, or the factorization problem mentioned above, the results were well 
constrained by the known measurement of ip' production, so the conclusion of the paper 
should hold. That is, the detection of rj' c and /i c ( 1 Pi) at HERA-P is feasible and promising. 
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TABLES 



TABLE I. Estimation of event numbers of /i c ( 1 Pi) and r]' c production at HERA-B 



State 






y/ 


bb 


inel. 


Cross section (/n) 


1076.1 nb 


153.7 nb 


84.2 nb 


12 nb 


13 mb 


Events rate (Hz) 


3311 


473 


259 


37 


40 M 


NP rod in 1Q 7 S 


3.3 x 10 10 


4.7 x 10 9 


2.6 x 10 9 


3.7 x 10 8 


4.0 x 10 14 


Final states 


7 J/V> 




IVc 




Vc + X 






Fraction 

(*+r 7 ( 7 )) 


1.2x 

io~ 5 


1.2x 
10~ 4 


1.2x 

io- 9 ~- 8 


1.2x 

10 -9~-8 


4.8x 

io- 8 


2.4x 

io- 7 




N prod in 1Q 7 S 

(*+r 7 ( 7 )) 


4.0 x 10 5 


5.6 x 10 5 


3 ~ 31 


3 ~ 31 


18 


89 




iV ofcs in 10 7 s 
(Assuming e = 10%) 


4.0 x 10 4 


5.6 x 10 4 


0.3 ~ 3 


0.3 ~ 3 


1.8 


8.9 





FIGURE CAPTIONS 

Figure 1. The leading order r( c production process at PN collision in both a s and v 2 . 

Figure 2. The generic diagrams of h c ( l Pi) production process at PN collision at leading 
order in v 2 ; (a) the Color-Octet process, (b) the Color-Singlet process. 

Figure 3. The generic diagrams of ifj' production process at PN collision; (a) and (c) the 
Color-Octet processes at v A and leading order in a s , (b) the leading order Color-Singlet 
process in both a s and v 2 . 
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